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and quantum fields satisfy. In this article we highlight some of the properties 
of the solutions to the equations and show to what extent they support the 
original ideas. 
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1. INTRODUCTION 

There is a tradition of condensed matter physicists and quantum field 
theorists (particle physicists) exploring ideas and techniques in common. 
The development of equilibrium renormalisation group methods in the '70s 
and onwards by both communities to describe phase transitions has been 
one of the major successes of that period. More recently, the adoption of 
simple causal bounds to constrain the non- equilibrium dynamics of continu- 
ous phase transitions has led to a renewal of this dialogue, reflected in the 
title Cosmological Experiments in Condensed Matter Systems of the review 
article by Zureka, the first extensive presentation of these ideas. 

Although subsequent work by several authors (including Zurek himselfiJi) 
has refined these original proposal^! they are, at heart, very simple: after a 
continuous transition the order parameter field (or fields) cannot adapt to 
a single uniform ground state value immediately. The reason is straightfor- 
ward. Although the adiabatic (equilibrium) correlation length £ e(J (i) diverges 
at the transition the true correlation length £(i) does not, since there is not 
enough time for it to do so. Causality imposes a maximum rate at which the 
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correlation length can grow and hence a maximum correlation length, £ say, 
at the onset of the transition. At the same time, causality imposes a horizon 
outside which the fields are uncorrelated. A consequence of a correlation 
length that is always finite is the creation of topological defects (monopoles, 
vortices, walls, etc.) that mark 'domain' boundaries. These defects then 
self-interact and annihilate, leading to larger and larger regions over which 
the field assumes one of the possible ground state values. 

If the simple relationship £ = 0(£def) between the correlation length £ 
and the separation length £^e/ of defects suggested by this picture is valid 
then the density of defects is bounded, and calculable, at their moment of 
formation. If their evolution is known, the density of defects at late times is 
similarly constrained. 

Since causality is as equally embodied in quantum field theory (QFT) as 
in condensed matter, these ideas had been posed independently by KibbleEm 
in the context of phase transitions of QFT in the very early universe. One 
of the great successes of QFT has been the unification of the electroweak 
forces through spontaneous symmetry breaking. There is every reason to 
believe that this and other symmetries were not always broken but that, in 
the very early and hot universe, they were restored. KibbleQ'u and othersa 
have observed that the same causal arguments put useful constraints on the 
density of defects at the time of their formation, which could have conse- 
quences today. Unfortunately, because of our lack of understanding of .the 
details of the early universe it is impossible to make reliable predictions^. 

It was Zurek who argued that the same causal bounds could be tested 
directly in condensed matter systems, whose phase transitions produce de- 
fects whose densities can be measured fairly readily. In the final section of 
this paper we review most of t h e experimental evidence for Zurek's bounds 
which, on several occasions! 1 1 H 1 2 H1j , are satisfied at a good qualitative level. 

Prior to that, in the next two sections we reiterate and rephrase the 
Zurek-Kibble bounds before providing an alternative scenario for the way 
defects appear after a transition. This is based on the fact that the simple 
defects that we shall consider have 'false' ground-state or vacuum at their 
cores where the field vanishes. Under suitable circumstances the separation 
length t^def is more sensibly derived by counting zeroes of the field as £<ie/ = 
0(£ 2ero )> where £ 2ero measures the separation of field zeroes. However, in 
principle, £ 2ero and £ are different correlation lengths, exploring different 
attributes of the fluctuation spectrum. 

We then apply these ideas in turn to condensed matter systems, using 
the empirical time-dependent Ginzburg Landau (TDGL) theory as specifying 
the dynamics, and to relativistic Quantum Field Theory (QFT). In the latter 
case, there is a further complication in that the system needs to decohere 
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(i.e. lose its quantum mechanical interference) before we can use classical 
probabilities to identify defects, and we devote Section 6 to that. 

Nonetheless, subject to certain conditions, we find that ^def is, indeed, 
0(£), as predicted, both for condensed matter and QFT. The Kibble-Zurek 
bounds are qualitatively valid, albeit for somewhat different reasons than 
simple causality. A byproduct is that the failure of 4 He experiments to give 
agreement0 may lie in an incorrect assumption of the decay rate of defects. 

We have not attempted to be exhaustive in this article and several 
important topics have not been included (e.g. temperature inhomogeneities) . 
The reader is referred to the references for more details. 



2. CAUSAL CONSTRAINTS ON THE PRODUCTION OF 

DEFECTS 

2.1. Condensed Matter 

The role of defects in our preliminary discussion was to both to infer and 
to predict the finiteness of the correlation length after the transition. Both 
condensed matter systems ( 3 He in particular^) and permit a prodi- 

gious variety of defects, and we shall report on some of them later. Rather 
than try to be generic we use the fact that most (but not all) experiments 
involve vortices and, with this in mind, we take an idealised superfluid as the 
prototypical condensed matter system in which to display these arguments. 
An essential property of superfluids is that their symmetries are global. For 
simplicity, we assume a single complex order parameter field </>(x, t), which 
permits the most simple vortices0. 

Most simply, the transition from normal fluid to superfluid is achieved 
through a pressure or temperature quench. The rate at which </>(x, t) can 
adjust, and the field become correlated, is the speed of (second) sound c(e), 
where e = (T — T c )/T c , which slows critically to zero as e — > 0. If we know 
the time dependence e(t) of (T — T c )/T c in the vicinity of the transition, the 
diameter of the 'sonic horizon' at time t is 

h(t) = 2 ( dt'c(t'), (1) 
J o 

where c(t) = c(e(t)), and we have assumed that the transition began at time 
t = 0, for which e(0) = 0. 

This immediately leads to a 'weak' causal constraint that, at time t 
after the transition has been completed, the correlation length £(i) satisfies 
the horizon bound 

m < h(t). (2) 



R.J. Rivers 



i.e. the field is uncorrelated outside its causal horizon. For long times @ 
becomes < 2c(ej), where ej is the final relative temperature. 

At his stage we need to be more explicit in that, for our prototypical 
case of a single complex field, there are two candidate correlation lengths, 
one for the field magnitude (e.g. superfluid density) and one for the phase 
correlation (e.g. superflow). Once the transition is complete it is the latter, 
associated with Goldstone modes, that determines the horizon. 

However, for the purpose of experiment the more relevant bound invoked 
by Zurek (and by Kibble in Ref.||) is the following 'strong' causal bound, for 
which this delineation is not necessary, in the first instance: Zurek argued 
that the correlation length freezes in before the transition, at t = —i, when its 
growth rate reaches its causal bound £ e <j(— t) = +c(— t), and only unfreezes 
at a comparable time +i after the transition. 

The major difference with the weak bound is that, before the transition, 
there is a single correlation length since the symmetry is unbroken. 
Further, as we shall see, the time t after the transition occurs early, just 
before the order parameter has reached its equilibrium value. At this time 
the real and imaginary parts of <f) are still effectively independent, and we 
still have only one correlation length. 

We reach the same conclusion from the viewpoint of the causal horizon. 
There is a time t, 

h(t) = (3) 

when the horizon is big enough to accommodate a defect. Before this time, 
< t < i , the adiabatic approximation has completely broken down, and 
the field is frozen in. Up to factors 0(1), this matches the value of i given 
in (||). It is argued that, at best, the system will emerge at time i with a 
correlation length £ = £ e q(t) that characterises identifiable domains. 

More specifically, in this and the other approaches, it is the longest 
wavelength modes, which control the field ordering, that freeze in first and 
unfreeze last. It must be stressed that these are rough estimates for the 
earliest time and the largest correlation length when the system permits an 
adiabatic description after the onset of the transition and we would not be 
surprised if, in practice, there was a difference by a factor of a few. 

Whatever, in a second assumption, that is independent of causality, 
Zurek further identified £ = S,eq(t) with £<fe/(*)> the defect (vortex) separation 
at the earliest time, t, at which they could be produced. To be specific, and 
simple, we assume that 

dt t q K J 



in the vicinity of t = 0, and adopt mean field critical indices for the equilib- 
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rium correlation length £ eq (t) and the speed of sound c(t), 



t 



'1/2 



c(t) = c 



t 



1/2 



(5) 



at this time, for suitable (zero-temperature) parameters £o and Co = £oAo- 
From ([|) it follows that 



and 



t ~ y/TQTQ = t Z 

1/4 



61 ? 



(6) 



(7) 



In practice, even for the fastest transitions, we have tq 3> to, whereby 
£,z S> £o- The variant of the causal argument that has the field frozen in with 
this correlation length at t = —iz suggests that « £z for — iz >t<tz- 
In particular £(0) « £z- 

From the second assumption we have that the initial vortex density n^ef 
=n de f(t) is 



ndef 



O 



pe \t q 



1/2 



(8) 



where / = 0(1) estimates the fraction of defects per 'domain'. Equivalently, 
the length of vortices in a box volume v is O^n^tv). 



2.2. Relativistic Field Theory 

Let us assume global symmetry breaking again. We note that if there 
had been no critical slowing down, as in relativistic QFT, and c(t) = c is 
constant, then the 'sonic' horizon h(t) of (Jxj) is replaced by the usual causal 
'light' horizon, to give 

e(i) < h(t) « 2ct (9) 

for the phase correlation length. This constraint was used to bound the 
production of monopolesl and cosmic strings! (vortices) in the early universe. 
On the other hand the strong causal bound (|3j) now gives 

Ik = i « (t q t 2 ) 1/3 (10) 

and 

^ = £ = £ e #i0~£o(^) • (ii) 
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In ( |TT|) is the common correlation length for early time, leading to the esti- 
mate fide/ = 0(1/£k) for the defect density at the time of their production. 

We stress that critical slowing down is not necessary for the existence of 
causal bounds. Its importance lies in the ameliorating effect it has on tem- 
perature inhomogeneities, whose effect is to diminish defect productionlll. In 
the early universe temperature inhomogeneity is slight, unlike in condensed 
matter physics experiments, where it is often substantial, and critical slowing 
down is crucial. 

This is all for flat space-time. 



2.3. Non-causal Mechanisms: The Ginzburg Regime 

In particular, what strikes us about the causal predictions is that they 
are as universal as the mean field approximation is valid. That, of course, 
is their main attraction. In the first instance the basic distance and time 
scales Co an d r o do not use any information about the magnitude of the 
order parameter (i.e. the strength of the interactions). Even the adoption 
of non-mean-field critical indices (necessary for 4 He, unnecessary for 3 He) 
hardly changes the picture. 

In stressing the primacy of causal horizons Zurek downplayed the im- 
portance of thermal fluctuations, whose characteristic scale is set by the 
Ginzburg temperature Tq through e(T(j), independent of the quench rate and 
dependent on the microscopic parameters of the theory. This is in contrast 
to the earlier suggestion (by Kibble! and others!) that thermal fluctuations 
in the Ginzburg regime might also lead to the production of vortices, again at 
early times. The reason why the Ginzburg regime might be important is the 
following: once we are below T c , the Ginzburg temperature Tq < T c signals 
the temperature above which there is a significant probability for thermal 
fluctuations between one degenerate groundstate and another on the scale 
of the correlation length at that temperature. That is, the thermal energy 
in such a fluctuation matches the energy required to pass over the hump of 
the unstable minimum. 

Whereas, above Tq one might anticipate a population of 'domains', 
fluctuating in and out of existence, at temperatures below Tq fluctuations 
from one minimum to the other become increasingly unlikely. Thus field 
configurations with non-trivial topology formed above Tq could stabilise. 
Since the Ginzburg regime does depend on the value of the order parameter, 
€.(Tq) can vary hugely from system to system. It had been suggested^ that 
we identify 

e eg (T G ) = , ^ (12) 
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with the scale £ at which stable domains begin to form. This is still invoked 
in cosmologyll. 

We shall see later why this is not the case. However, even if thermal 
fluctuations are not the mechanism for the formation of larger domains, they 
cannot be ignored completely. At the least they are relevant to the formation 
of small 'domains', and to wiggles in the boundaries of larger domains where 
defects are to be found. As such, they can make the definition of defect 
density scale invariant. 



3. AN ALTERNATIVE APPROACH: INSTABILITIES AND 

(LINE) ZEROES 

In the several years since these simple bounds were first proposed we 
have acquired a much better understanding of the way in which transitions 
occur. These does not mean that these bounds have lost their relevance, but 
that they need to be qualified. 

To be more quantitative we assume that our single complex scalar field 
(p in three spatial dimensions has a wine-bottle potential. The transition 
is continuous. That is, we assume that the qualitative dynamics are condi- 
tioned by the field's equilibrium free energy, of the form 

F(T) = J d 3 x (JV<A| 2 + e(T)H 2 + AH 4 ) . (13) 

where we have scaled the field to be as generic as possible. 

On rescaling eft could represent a relativistic quantum field, with free 
energy 

F(T) = J d 3 x (jV</>| 2 + m 2 (T)|0| 2 + AH 4 ) . (14) 

where m 2 (T) = M 2 e(T) at one-loop level. On the other hand, on scaling, F 
could be the Ginzburg-Landau free energy 

F(T) = J d 3 x (|^|V0| 2 + a(T)\<p\ 2 + /?H 4 ) (15) 

in which the chemical potential a{T) = aoe(T) vanishes at the critical tem- 
perature T c . With minor qualifications (|l5|) is the free energy adopted for 
simplified 3 He in Refjl^ and Bose-Einstein condensation^ (BEC). It is not 
a reliable representation of 4 He, but may still be useful in this case. 
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3.1. Time-scales: The Spinodal Time 

Any dynamical equations for the onset of a continuous transition will 
embody causality, by definition. However, the transition cannot be said to 
have happened before the order parameter has achieved its equilibrium value 
\(j)\ 2 = 1/A, (or M 2 /X or ao/(3, depending on how we rescale (|T^)). If (...)( 
denotes ensemble averaging at time t then a lower bound on the first time 
from which we can start counting defects is t = t sp , for which (\4>\ 2 )t = 1/A. 

Ensemble averaging is taken with respect to the relative probability pt[(p} 
that the order parameter field 4> takes the value ^>(x) at time t. The way in 
which (\4>\ 2 )t builds up to its final value is by the growth of the amplitudes 
of the unstable long-wavelength modes, which are unstable because of the 
upturned parabolic free energy at initial times. These long wavelength modes 
order the field on increasingly larger scales. The time t sp is, crudely, the time 
for these modes to roll from the top of the hill to the groundstates at the 
bottom, and we have termed it the spinodal time. This is a very different 
picture from that of a field freezing in by, or before, the transition that was 
proposed in Section 2. 

A priori, t sp is not related to the causal ix or iz, but it is not difficult 
to see why they might be comparable. Unstable modes grow exponentially 
fast. As long as dimensional analysis makes ix or iz the natural unit in 
which to measure time, any exponentially growing term will achieve values 
that are not exponentially large at times t = 0(1) in these units. We thus 
anticipate only a logarithmic sensitivity to the microscopic parameters of 
the theory (and the quench rate). 

3.2. Quantum Field Theory: The Decoherence Time 

Although the spinodal time is equally important for QFT, there are 
additional problems in that the evaluation of ensemble averages (...)f is not 
enough. In particular, in QFT we need to consider the whole density matrix 
(4> + \p(t)\(j)^) rather than just the diagonal elements pt[<fi] = ((f)\p(t)\(p) that 
determine the probability functional. 

The reason is that probabilities are only useful when there is no, or 
little, quantum interference between adjacent configurations (cf. the two-slit 
problem) . Such interference is measured by the magnitude of the off-diagonal 
density matrix elements. 

Before we can think of identifying particular field configurations like vor- 
tices in QFT we must have that the reduced density matrix is approximately 
diagonal on coarse-graining. 

By coarse-graining we now mean the separation of the whole into the 
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'system', the long wavelength modes of which establish the field ordering, 
and its 'environment' whose degrees of freedom are integrated over to pro- 
duce the reduced matrix. The environment comprises the other fields with 
which our scalar is interacting, together with the short wavelength modes of 
the scalar field itself EE. This coarse-graining is more sophisticated than 
the effective short-distance cutoff that is essentially all that is needed in 
condensed matter. 

The reduced density matrix is non-unitary (dissipative) and the effect 
of this environment is to introduce a new time scale, the 'decoherence' time 
such that quantum interference can be ignored at times later than it. 
Only for t > tn does it make sense to think of classical defects. However, 
dissipation will also grow exponentially initially since it is driven by the same 
long wavelength unstable modes. Thus, although to depends on interaction 
strengths as well as the quench rate, we anticipate that it is insensitive to 
them, for the same reasons that i is insensitive. For the moment we assume 
that tjj <t sp to keep the discussion simple. 



As we noted earlier, this dynamic picture, in which modes build up their 
amplitudes from their initial small Boltzmann values when the field is sitting 
at the top of the hill, by rolling down it, or spreading over it, seems totally 
at variance with that given in Section 2. There, it was assumed that the 
field freezes in at, or even prior to the transition, with a correlation length 
£ (£,K or £z) that is already huge, comparatively. 

In fact, even at early times when we can get by with a single correlation 
length for the 0-field, we are talking about two different lengths. The length 
£ is the usual correlation length, obtained from the long-distance behaviour 
of the correlation function at time t, 



(r = |x|) in which the P(k,t) is the power spectrum of the fluctuations. To 
obtain £ we take r — ► oo, and look for behaviour of the form exp— (r/£) 7 . 
That is, £ is determined by the nearest singularity of P(k, t) in the fc-plane. 

This is not the length that characterises the separation of defects, which 
are formed from the phase separation as the field falls from the hill in different 
directions at different spatial points. Vortices form because these phases do 
not match up smoothly and are frustrated by zeroes in the field. 

More specifically, at long times after the transition, we expect to find 
widely separated classical vortices. The classical vortices of this simple the- 
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(16) 
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ory, solutions to the equation 8F/S(j> = are tubes of false groundstate, 
with line-zeroes eft = at their cores, with width 0(£ eq (T)). If we write 
cj) = (0 1 + 4>2)/V^ these line-zeroes are the intersections of the surfaces 
4>\ = 4>2 = 0. From the earliest numerical simulations onwardsll@ simple 
vortices (and defects in lower dimension) have been identified by counting 
zeroes of the fields. 

That is, the empirical measure of defect separation is £ 2ero , the typical 
separation of line-zeroes. If, for the sake of argument, we assume Gaussian 
field fluctuatiDJis, then we shall see that, at their time of formation, ^ zero (t) 
is given byc3'E3 

Primes denote differentiation with respect to r. Empirically!!, ( |i~7|) has been 
known to be valid until the spinodal time, even though the field has ceased 
to be Gaussian by then. 

This suggests that a more realistic length-scale for calculating defect 
densities at the onset of defect production is £ 2ero = £, ze ro(tsp)- In this case 
there is no immediate reason to relate £ 2ero to £ since, unlike £ e(3 (i), £, ze ro(t) 
is given entirely by the s/iori-distance behaviour of G(r, t) or equivalently, 
by the dominant wavenumbers in P(k,t). 

It is not difficult to extend the picture of vortices as line zeroes to other 
defects like global monopoles and domain walls in three dimensions, or kinks 
in one dimension. No novelty arises, and we shall not do so. 



3.4. When are Line Zeroes Vortices? 

We have rather assumed that, even if we cannot identify £, zero with £ eq (t) 
(i z or i K ), we can identify £ zero with £ de f. 

This cannot be the case exactly in that, whereas all such defects can 
be identified by their line-zeroes (or zeroes), not all line zeroes (or zeroes) 
are candidate defects, since zeroes occur on all scales. A starting-point for 
counting vortices is to count line zeroes of an appropriately coarse-grained 
field, in whichstructure on a scale smaller than £q> the classical vortex size, 
is not presented. Lattice-based numerical simulations do this automatically 
(but see Ref.^). For the moment, we put in a cutoff I = O(£o) by hand. 
We note that the inclusion of a cut-off does not affect the long-distance 
correlation £. 

Suppose that, at time t sp , the typical separation of line-zeroes for this 
coarse-graining is ^ ze ro,i(t sp )- Before line zeroes can be identified with clas- 
sical (global) vortex cores, and £ Z ero,i(tsp) with Cdef(tsp), the following con- 
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ditions are necessary. 

• Although the correlation length is insensitive to a short-distance 
cutoff Z, this is not the case for £, zero ,i{t)- Thermal fluctuations will give 
structure on small scales which will lead to ambiguity in its definition 
and in the subsequent density. Only when d£ zer0j i/dl is small in com- 
parison to £zero,i/l at Z = £q wm the line-zeroes have the small-scale 
non-fractal nature of classical defects, although defects may behave like 
random walks on larger scales. Among other things, this will depend 
on the Ginzburg regime. 

• As a final, related, check, the energy in field gradients should be com- 
mensurate with the energy in classical vortices with the same density 
as that of line zeroes. 

In fact, most (but not all0) numerical lattice simulations cannot dis- 
tinguish between line-zeroes and classical vortices. 

We shall see that, roughly, fluctuations separate into thermal fluctu- 
ations controlled by the current temperature T(t), which determine the 
small-scale structure of line-zeroes, which compete with large-amplitude 
long- wavelength fluctuations whose role is to provide large-scale order. It 
is these latter with which classical defects are associated. 



3.5. Counting Zeroes 

Suppose, at some time, that the field has line zeroes R„,(s), where n = 
1,2, .. labels the zero, and s measures the length alone it. As a result the 
topological line density of zeroes p(r) can be definedc3u3 by 

p {x ) = J2jds^5^-K n (s)]. (18) 

In (jilj) ds is the incremental length along the line of zeroes R n (s) and 
dR n /ds is a unit vector pointing in the direction which corresponds to pos- 
itive winding number. 

It follows that, in terms of the zeroes of the field </>(x), pi(x) can be 
written as 

Pi(x) = (5 2 [0(x)]e ijfc aj^i(x)a fc (/>2(x), (19) 

where <5 2 [(/>(x)] = 5[</>i(x)]<5[(fo(x)], where (j) = {<j)\ + i(j)2)/V2. The coefficient 
of the 5-function in ( |l9| ) is the Jacobian of the more complicated transfor- 
mation from line zeroes to field zeroes. What we want is not this, but the 
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total line density p(x), 

Pi{x) = 5 2 [0(x)]|e Jjfe a j 1 (x)9 fc( / )2 (x)|. (20) 

Let us suppose that, at time t, the probability density that the field 
takes functional form 0(x) is Pt[<^>]. The ensemble average of K[cp] at time t 
is 

{K[(t>]) t = J VfaVfoPtW *[</>]■ (21) 

The vanishing field expectation value and the independence of the field and 
its derivatives 

(M*))t = = ^„(x)a^ 6 (x)) tj (22) 

imply (pj(x))t = i.e. an equal likelihood of a string line-zero or an anti- 
string line-zero passing through an infinitesimal area. However, 

(t) = (pj(x)) t > (23) 

and measures the total line-zero density in the direction i, without regard 
to string orientation. The isotropy of the initial state guarantees that n is 
independent of the direction i. 

To get a basic idea as to what this means let us assume that the field fluc- 
tuations are Gaussian. This has been the starting assumption for networks 
of cosmic stringsll'B. Then everything is given in terms of the two-field 
correlator (16) at time t, 

(Mx)M0))t = 5 ab G(r,t). (24) 
We coarse-grain the field by putting in a cutoff I = O(£o) by hand, as 

Gi(r,t) = J fke ik - x P(k,t)e- k2l \ (25) 

where we assume I = O(£o) or 0(^ eq (Tf ina i)), as appropriate. 

As we indicated, 
termined completely^ 



ier, in this Gaussian approximation n zero j(t) is de- 
by the short- distance behaviour of G(r, t) as 



1 _ -l Gf(0,t) 
nzer ° At) ~ 2^ zero ^t) 2 ~ 2^ GKM' (26) 



where we have used (|26|) to define £ Z ero(t). That is, n zero {t) is determined 
by the ratio of the fourth to second moments of P(k,t). 

For non-Gaussian fields the situation is much more complicated. How- 
ever, as long as there is a dominant wavenumber ko(t) in P(k;t) this sets 
a length scale £ ~ /cq(^) _1 that characterises vortex separation. As we said 
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earlier, (|26| ) can be approximately valid until the spinodal time, and we shall 
assume that to be the case. 

We note that a line-zero is invariant under gauge transformations <j) — ► 
4>e ia . Whereas a simple spatial cutoff breaks gauge-invar iance, we have 
learned from lattice gauge theory how to coarse-grain a gauge field. However, 
for a local theory, taking 4>i and (j>2 independently Gaussian is a gauge- 
dependent statement. We shall not pursue this further. 

We conclude with a brief discussion of topological charge through a 
surface, for which the topological density (|H]), rather than the total den- 
sity (p0|), is appropriate, again in the Gaussian approximation. Consider 
a circular path in the bulk material (in the 1-2 plane), circumference C, 
the boundary of a surface S. For given field configurations ^ a (x) the phase 
change 9c along the path can be expressed as the surface integral 

9 C = 2vr f d 2 x p 3 (x), (27) 
Jxes 



Again we quench from an initiaLstate satisfying ( p2[ ) (i.e. no rotation). 
It is not difficult to showBll that, if f(r,t) = G(r,t)/G{0,t), then 

^-a^r**- <28) 

The linear dependence on C is purely a result of Gaussian fluctuations. We 
note that (A#c) 2 requires more than the very short distance behaviour of 
fi. Nonetheless, if P(k,t) is strongly peaked at k = ko, then £ s = 0(kQ 1 ), 
as is ^zero- Further, if we removed all material except for a strip from the 
neighbourhood of the contour C we would still have the same result. This 
supports the assertion by Zurek that the correlation length for phase varia- 
tion in bulk fluid is also appropriate for annular flow. 

That (A6c) 2 varies in time is a consequence of defects migrating through 
the boundary from the inside to the outside and vice- versa, and annihilating. 
For a finite system, where such migration is impossible and topological charge 
is conserved, we would think of adopting ( pgj ) at the moment that defects 
appeared. 



4. WHEN DO THE APPROACHES AGREE? CONDENSED 

MATTER 

Let us summarise the situation for our idealised condensed matter, as 
far as the 'strong' causal bounds of Zurek and Kibble are concerned.. The 
time and length scales of Zurek's analysis are tz of (|6|) and £z = ieq{tz) 
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of (|7])- In our analysis of zeroes the relevant time and length scales are t sp 
and £, Z ero = £,zero(t S p)- Above and beyond that, thermal fluctuations may 
complicate the identity of ^ ze ro and Cdef- 

We have already suggested that tz and t sp may be comparable, quali- 
tatively, but the situation for length scales is less simple, not least because 
of thermal fluctuations. To make further progress requires a concrete model 
in which correlation functions can be calculated explicitly. 



4.1. The Time-Dependent Ginzburg-Landau (TDGL) Equation 

Not surprisingly, because the Time-Dependent Ginzburg-Landau (TDGL) 
equation lends itself so directly tonumerical analysis, it has been the basis 
of many papers, largely by ZurekErH and co-workers (but not exclusively e.g. 
see Refs jl^pO| ). By definition it encodes causality. In Ref.|i~9| a stripped- 
down global U(l) model was used to mimic the Helsinki 3 He experimentta. 
A more realistic TDGL model for 3 i^e^that takes the full order-parameters 
into account has been used elsewhere^. The TDGL equation is also useful 
in BEC§. 

In the last resort, only numerical methods permit us to get long-time 
solutions, but in this article we pick out some properties of simplified analytic 
solutions that we think show the relevance of the causal bounds. 

We assume that, for the condensed matter systems of interest to us, the 
dynamics of the transition can be derived from the explicitly time- dependent 
Landau-Ginzburg free energy (£q = 1) 

F(t) = J d 3 x Q(V^) 2 + ^(t)^ + I A (^) 2 ). (29) 

in which we substitute T(t) for T directly in (|l~3|), where (j) = ((/)% + ifa) / V% 
(a = 1,2) is the complex order-parameter field, 

For a simple dissipative system the time-dependent Landau-Ginzburg 
(TDLG) equation is 

r^r = -Wa +Va ' (30) 

where rj a is Gaussian thermal noise, satisfying 

{ Va (x, t)rib(y, t')) = 2d ab T(t)rS{ X - y)S(t - t'). (31) 

The basic timescale to, the relaxation time of the long wavelength 
modes, is T _1 in our units. 

In terms of the dimensionful constants of ([D]) the fundamental length 
scale £o is £o = fr 2 /2mao whereby Ceq(T) = £q(T/T c — l) -1 / 2 as before. With 
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ro = l/r«o it follows that the equilibrium correlation length £ eq (t) and the 
speed of sound behave when t vanishes as in (||), 



~ 1/2 -m ^ 

, c(t) = — 

^0 



1/2 



(32) 



from which the causal bounds ([]) and (]?]) follow. 



4.2. Preezing-in Happens (and is Irrelevant) 

It is relatively simple to show that the TDGL equation (^) embod- 
ies Zurek's causality bound. At early times \t\ < iz the effective potential 
V((p, T) is still roughly quadratic and the self-interaction term can be ne- 
glected (A = 0). 

In space, time and temperature units in which £o = ro = ks = 1, Eq.pOj 
then becomes 

a (x,t) = -[-V 2 + e(t)]^ a (x,t) + f/ a (x,t). (33) 

where fj is the renormalised noise. The solution of the 'free'-field linear equa- 
tion is straightforward, giving a Gaussian equal-time correlation function 

G{r,t) = J fkj*- T P(k,t). (34) 

in which the power spectrum P(k, t) has a representation in terms of the 
Schwinger proper-time r as 

P(k, t) = [°° dr f{t - r/2) e~ Tfc2 e~ £ ds e(t_s/2) , (35) 
Jo 

where T is the renormalised temperature. In turn, this givesEHIl 

G(r, t) = J™ dr f(t - r/2) e " Jo ds < t ~ s ^ . (36) 

For e(t) of Eq.|] a saddle-point calculation shows that, at time t = 0, 
when the transition begins, 

G ( r7 0)«^L e -^/^) 4/3 , (37) 

on rescaling, where a = 0(1), confirming Zurek's result of a frozen field 
correlated over many cold defect thicknesses. 
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It is immediately apparent that £z) is n °t a measure of the separation 
of line zeroes. Putting in the momentum cutoff A; -1 > I = Z£o = O(£o) °f 
Eqj3q by hand corresponds to damping the singularity in G(r, t) at r = 
aso 

making Gi(0,t) finite. We stress that, for t « 0, the correlation length £ 
remains 0(£), independent of Z. 
At i = both Gj(0,t) and 

- G "( o •^ 2 '^^p|^ e - /;J "' , '-* /2, ■ < 39 > 

are dominated by the short wavelength fluctuations at small r. Even though 
the field is correlated over a distance £z 3> Z the density of line zeroes 
n zero = 0{l~ 2 ) depends entirely on the scale at which we look. Equivalently, 
Czero(O) = 0(1). In no way would we wish to identify these line zeroes with 
prototype vortices and £ 2ero (0) has nothing to do with £z at this time. 



4.3. Why t sp « t z and £ zero « £ def r* £ z (in General) 

As the system evolves away from the transition time, the free equation 



Eq.33 ceases to be strictly valid but, to a first approximation, the back- 
reaction does not set in until the field has sampled the groundstates. We 
can keep the linear approximation, with its single correlation length, until 
t ~ t sp . For the linear quench of (|4j) we findE 

irnt-r/2) [4j(r + i2)]3/2 ■ (40) 

However, as time passes the peak of the exponential grows and n zero be- 
comes increasingly insensitive to Z. How much time we have depends on the 
magnitude of A -1 , since once G(0, t) has reached this value it stops growing. 
Provided that the peak at r = 2t dominates over the thermal fluctuations 
at r ~ then 

_ „ /-oo „-(r-2t) 2 /4p„ 

Since G(0, t) = 0(exp((t/iz) 2 ) at early times the backreaction is imple- 
mented extremely rapidly, justifying the free-field approximation in the ex- 
pression for Gi(0,t) above. 
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Let us assume that the thermal fluctuations can be ignored by time 
t = iz- In our units, in which XT C = y/1 — Tg/T c , the condition that 
{\(j)\ 2 )t = 1/A at t = t B p gives 

hz ~ I L ( t q/ t ° 

tz~2V \(1-T G /T c f 

Although greater than unity, the ratio is O(l), extremely insensitive to the 
parameters of the model, and to the quench rate. For systems and experi- 
ments as widely different as the ones for 3 He and 4 ffe that we shall discuss 
later, we have t sp /tz ~ 3. 

If we continue to assume that the thermal fluctuations at r ~ can be 
ignored, then we can use (|3S|) and (j39|) to calculate the density of line zeroes. 
The result is 

t 1 

in accord with (H) of Zurek, where f 2 8irt sp /tz = O(10 2 ). Equivalently, 
the line-zero separation is 

Izero « e 2 ero(t S p) = 2(t sp /i Z ) 1/2 e = 0(fz), (44) 

even less sensitive to the parameters of the model and the quench rate. 

By assumption, n zero is independent of the cutoff I (I — O(l) in our 
units). The line-zeroes are on the verge of being classical vortices since, 
once (|0| 2 )t = 1/A, the Gaussian field energy, largely in field gradients, is 

F*>(J d 3 x i(V<^) 2 ) = -VG"(0, t sp ), (45) 

where V is the spatial volume. This matches the energy 

E « Vn def (t)(2<irG(0, t sp )) = -VG"(0, t sp ) (46) 

possessed by a network of classical global strings with density n zero , in the 
same approximation of cutting off their logarithmic tails. 





4.4. Why Back-reaction may not Matter (Much) 

The full TDGL equation is 

<^ a (x, t) = -[-V 2 + e(t) + A|0(x, t)| 2 ]</) a (x, t) + f/ a (x, t). (47) 
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In order to retain some analytic understanding of the way that the den- 
sity is such an ideal quantity to make predictions for, we adopt the approxi- 
mation of preserving Gaussian fluctuations by linearising the self-interaction 

as 

0«(x,i) = -[-V 2 + e e// (t)]<Aa(x,t)+f/ a (x,t), (48) 

where e e // contains a (self-consistent) term 0(X(\(j)\ 2 }). Additive renormali- 
sation is necessary, so that e e ff « e, as given earlier, for t < to- 

Self-consistent linearisation is the only usable approximation in non- 
equilibrium QFTE3, but is not strictly necessary here, since numerical sim- 
ulations that identify line zeroes of the field can be made that use the full 
self-interaction@'H. However, to date none address the questions we are pos- 
ing here exactly, and until then there is virtue in analytic approximations 
provided they are not taken too seriously. Our discussion complements that 
of Ref.^8| in which a self-consistently linearised TDGL theory is also exam- 
ined. 

The solution for Gi(r,t) is a straightforward generalisation of (36), 

making Gi(0,t) finite. 

Assuming a single zero of e e ff(t) at t = 0, at r = the exponential in 
the integrand peaks at r = f = 2t. Expanding about f to quadratic order 
gives 

„ ft , I n| «o -(r-2i) 2 |e'(*o)]/4 

Gl (0,t) « ^/o^/MI jT drT{t-r,2) — - p)f/2 . (50) 

The effect of the back-reaction is to stop the growth of Gi(0,t) — G;(0,0) = 
(l^| 2 )t ~~ (l^| 2 )o at its symmetry-broken value A -1 in our dimensionless units, 
thereby preserving Goldstone's theorem by requiring \e e ff(u)\ — ► oo. 

What is remarkable in this approximation is that the density of line ze- 
roes uses no property of the self- mass contribution to e e ff(t), self-consistent 
or otherwise. With 

« 2^e 2 Jo du l^//( u )l / drT{t-rl2)—— l ^7^(52) 







[4vr(r + P)fl 2 



all prefactors in Ti zero cslucgI, to giv> 



nZer0{ ) " 4- J- j^fit - r/2) e -(-(-o))74*l ' ^ 
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thereby justifying our free-field approximation for t < t sp . It is for this 
reason that simple dimensional analysis (the basis of the causal bounds) is 
so successful. 

For t > t sp the equation for n zero (t) is not so simple since the estimate 
above, based on a single isolated zero of e e //(i), breaks down because of the 
approximate vanishing of e e ff(t) f° r t > t sp . A more careful analysis shows 
that Gi(0,t) can be written as 

^-r iffin^ ^' (54) 

where G(t, t) has the same peak as before at r = 2t, in the vicinity of which 

G( T ,t) = e 2 io* **l«e//(«)l e -(r-2t)V4^ > (55) 

but G(r,t) = 1 for r < 2(i - i sp ). Thus, for t q > r , G;(0,t) can be 
approximately separated as 

G l (0,t)^GY v (t) + G IR (t), (56) 

where 

roo 

GY v {t) = / drf(t-r/2)/[47r(r + P)] 3/2 (57) 
Jo 

« T(t) / dr/[47r(r + P)] 3/2 (58) 
Jo 

describes the scale- dependent short wavelength thermal noise, proportional 
to temperature, and 

1 f°° 

GlR{t) = (8^ J jrf{t-T/2)G{T,t) (59) 



describes the scale- independent, temperature independent, long wavelength 
fluctuations. The integral /f^, drG(r, t) is naturally time-dependent, largely 
cancelling its prefactor so as to keep (\<fi\ 2 )t constant. 

A similar decomposition G///(0,t) Gf/f v (t) + G// IR (t) can be per- 
formed as 

roo 

Grr¥ v (t) = 2vr / dr T(t - t/2)/[4tt(t + P)] 5 / 2 (61) 
Jo 

roo 

« 2vrT(t) / dr/[47r(r + P)] 5 / 2 . (62) 
Jo 
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and 

G" IR {t) = —-^ J JrT(t-T/2)G{T,t) (63) 



(87rf) 5 /2 

47rf c 
(4vrr(t)) 5 / 2 



-oo 



drG(r,t). (64) 



In particular, Gff IR (t)/G IR (t) = O^ 1 ). 

Firstly, suppose that, for t > t sp , G IR (t) > Gf v (t) and Gu IR {t) > 
Gf/f v (t), as would be the case for a temperature quench T(t) — * 0. Then, 
with little thermal noise, we have widely separated line zeroes, with den- 
sity n zero (t) bs -G// IR (t)/27rG IR ( y t). With dn zero/ 91 small in comparison to 
n ze ro/ 1 at I = £o we identify such essentially non- fractal line-zeroes with pro- 
totype vortices, and n zero with ridef- Of course, we require non-Gaussianity 
to create true classical energy profiles. Nonetheless, the Halperin-Mazenko 
result may be well approximated for a while even when the fluctuations are 
no longer GaussianE3. This is supported by the observation that, once the 
line zeroes have straightened on small scales at t > t sp , the Gaussian field 
energy continues to match that possessed by a network of classical global 
strings with density n zero . 



4.5. The Role of the Ginzburg Regime 

We have seen that the UV thermal fluctuations at time t come from 
the small r part of the integration and are, approximately, proportional to 
T(t). It is the thermal fluctuations that give rise to small-scale structure 
on the defects and prevent us from giving a scale-independent value to their 
density. Thus, if we reduce the temperature to absolute zero, the thermal 
fluctuations vanish and the density becomes well defined. 

Our first observation is that, in practice, this never happens. In pressure 
quenches it is the critical temperature that changes, at almost constant 
temperature. Even in cooling, the final temperature is usually a substantial 
fraction of the critical temperature. 

For the sake of simplicity, let us keep Tf = 0(T C ) < T c and thereby 
take T = T c in Gi(0,t) above. The necessary time-independence of G IR (t) 
for t > t sp is achieved by taking e e ff(u) = O^" 1 ). In consequence, as t 
increases beyond t sp the relative magnitude of the UV and IR contributions 
to Gi(0, t) remains approximately constant. Further, since for t = t sp , 

e 2/;^w«)i e -(At)7*i^i, (es) 

this ratio is the ratio at t = t sr) . 
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Nonetheless, as long as the thermal fluctuations are insignificant at t = 
t sp the density of line zeroes will remain largely independent of scale. This 
follows if G" IR (t sp ) > G< y (t sp ), since GV/;(0,t) becomes scale-independent 
later than G[(0,t). In l£ we showed that this is true provided 

(t q /t )(1-T g /T c ) <Ctt\ (66) 

where C = 0(1) and Tq is the Ginzburg temperature. That is, for slow 
quenches or quenches in which thermal fluctuations remain large, there are 
no identifiable defects at the spinodal time, since the line zeroes are highly 
fractal on small scales and the Zurek analysis breaks down. 

In our self-consistent linearisation the situation never improves, al- 
though this is extreme. Of course, if we view line-zeroes through a lattice, 
they will be seen. Choose a different lattice, and we will see more, or less. 
Numerical simulations that identify vortices with line zeroes (and nothing 
more) on a fixed lattice are suspect until thermal fluctuations become unim- 
portant. Certainly, for uniform quenches, defects are not formed at the end 
of the Ginzburg regime by thermal fluctuations freezing in, as originally 
suggested. 

A simple way to see if defects can be produced by thermal fluctuations 
is to keep the system in the Ginzburg regime for a long period At 3> t sp and 
then drop out of it. Numerical simulations have been performed by Zurek 
and collaborators^ that test this possibility. Empirically, the end result 
is not very different from the case of the linear quench that we have been 
discussing hitherto. The thermal fluctuations of the Ginzburg regime are 
not the source of defect production. 

To get an analytic idea as to why this is the case we return to the 
linear approximation of the previous section. Even with more complicated 
temperature profiles the separation into long and short wavelength modes 
will still occur. 



In particular, we expect the end-point behaviour (57) and ( pT| ) for the 
scale-dependent short wavelength fluctuations to be unchanged. As a result, 
Gy v (t) and G/t¥ v (t) are largely insensitive of the past history of the system, 
even if that history involved a long time At in the Ginzburg regime. 

Suppose, as in Ref J4C|, we begin at temperature To, then drop to tem- 
perature Tq for a period At, and then drop to temperature T = 0. In 
comparison to spending the whole time at Tq, the strength of the short 
wavelength fluctuations in G//Y v (t) is 

K^(I-t))- 

for t >> At, with comparable behaviour for Gj'(t). That is, if we end up 
at low temperature, our intermediate history is largely unimportant for the 
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small-scale fluctuations, however that temperature is arrived at. In partic- 
ular, short-distance scale-dependent fluctuations will be negligible at long 
times, whatever. 

Similarly, provided the temperature is monotonically decreasing (since 
otherwise the effect on G\ R (i) and Gll\ (t) is dramatic) we do not expect a 
strong effect of past history on the IR contributions beyond a delay in the 
position r of the peak in r by approximately 2Ar. Since f(t) depends on 
where e e ff(t — t/2) vanishes, we must have f(t) = 0(t), largely insensitive 
to intermediate conditions, once t is large enough. On neglecting the (now 
negligible) UV terms, we have n zero determined approximately by f(t) alone 
(giving the usual n zero oc t^ 1 behaviour V To a good approximation this is 
what is seen in the numerical simulations^. To quantify the slight differences 
seen in the numerical simulations arising from different temperature profiles 
requires further work. 

There is no way that defects can be produced by attempting to incubate 
them in the Ginzburg regime. 

5. WHEN DO THE APPROACHES AGREE? RELATIVISTIC 

QFT 

For QFT the situation is rather different. In the previous section, in- 
stead of working with the TDLG eauation, we could have worked with the 
equivalent Fokker-Planck equations for the probability pt[&] that, at time 
t > 0, the measurement of (ft will give the function ^(x). When solving the 
dynamical equations for a hot quantum field it is convenient to work with 
probabilities from the start. 

5.1. The Closed Time-path 

Take t = to as our starting time for the evolution of the complex field 
(ft = ((f> i +i(f>2)/V2. Suppose that, at this time, the system is in a pure state, 
in which the measurement of (ft would give <I>o(x). That is:- 

4>(t = 0,x)|$ ,t = to) = $o|$o,t = t ). (68) 

The probability pt[3>] that, at time tf > to, the measurement of (f> will 
give the value $ is pt[$] = l^ol 2 ! where is the state- functional with the 
specified initial condition. As a path integral 



* = / V(t>exv{iSt[(f>}\, (69) 
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where St[4>] is the (time-dependent) action that describes how the field is 
driven by the environment and spatial and field labels have been suppressed 
(e.g. Vcf) = V(f)iD(f)2)- Specifically, for t > to the action for the field is taken 
to be 

SM = J dx^-d^d^a - \e(t)M 2 ^ a - ^A(^) 2 ). (70) 

where e(t) is as before. Henceforth we return to natural units in which £o = 
r = M- 1 = 1. 

It follows that pt[&] can be written in the closed time-path form 

PM = T ^ V^V4>- exp{*(s t [0 + ] - S t [<f>-])\, (71) 

where = Dcpf'Dcjyf' . Instead of separately integrating c/> ± along the time 
paths t < t < tf, the integral can be interpreted as time-ordering of a field 
4> along the closed path C + © C_ where 4> = 4> + on C+ anci 4> = on C- • 
When we extend the contour from tf to t = oo either or 0~ is an equally 
good candidate for the physical field, but we choose <p + . 

The choice of a pure state at time t = to is too simple to be of any use. 
Let us assume that $ is Boltzmann distributed at time t = to at an effective 
temperature of Tq = (5q 1 according to the Hamiltonian H [$] corresponding 
to the action S [(/>], in which </> is taken to be periodic in imaginary time with 
period Pq. We now have the explicit form for pt[G>], 

p t [$] = f V4>e iS ^ S^+itf) - $], (72) 

written as the time ordering of a single field along the contour C = C+ffi(7_ffi 
C3, extended to include a third imaginary leg, where takes the values <p + , 
4>~ and (ps on C+, C_ and C3 respectively, for which Sc is 5'[0 + ], S[4>~] and 
So [(^3]. It is sufficient for our purposes to take the Boltzmann distribution 
on C3 to be the free- field distribution with mass parameter m 2 = M 2 eo, 
corresponding to a temperature quench from an initial state at temperature 
T > T c , where (T /T c - 1) = e . 

As with condensed matter, it is not necessary to calculate pt[$] directly 
since G a fe(|x - x'|;i) = ($ a (x)$ 6 (x')) t is given by 

G a6 (|x-x'|;t) = ^ a (x,t)^(x',t)), (73) 

the equal-time thermal Wightman function with the given thermal boundary 
conditions. 



*0 
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Fig. 1. The closed timepath contour C+ © C_ with a third imaginary leg 



5.2. The (Irrelevant) Freeze-in 

Unlike for condensed matter we shall not postulate dissipative equations 
(with corresponding noise), but just calculate the evolution of the quantum 
field under its time-dependent Hamiltonian derived from (|29|). Insofar as 
the results look similar it is for totally different reasons. [For simplicity we 
continue to work in flat space-time.] 

As for the condensed matter case, the interval — tn < t < ijc occurs in 
the linear regime, when the self-interactions are unimportant. The relevant 
equation for constructing the correlation functions of this one-field system 
is now the second-order equation 



5F 



dt 2 



(74) 



for F of Eq.13. This is solvable in terms of the mode functions X^if), 
identical for a = 1,2, satisfying 



dt 2 



+ k 2 + e(t) 



(75) 



subject to Xk (t) = e^' iU)int at large negative times, with incident frequency 
uoin = Vk 2 + e . 



The diagonal correlation function G(r, t) of Eq.34 is given as the equal- 
time propagator 



G(r,t) 



k^xUt) x - k (t)c{k) 



(76) 



where C{k) = coth(ui n (k) /2Tq) /2u>i n (k) encodes the initial conditions. 
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An exact solution can be givene3 in terms of Airy functions. Dimensional 
analysis shows that, on ignoring the k-dependence of C(k), appropriate for 



large r (or small k), £k = £,eq(ti<) of Eq|ll| again sets the scale of the equal- 
time correlation function. Specifically, 

e(,,0)«/**5&#F W , (77) 

where F(0) = 1 and F(k) ~ for large k. 

Kibble's prediction for the large scale at which the field freezes in is 
correct but, as with condensed matter, has nothing to do with the formation 
of defects around t = 0. At t = there are, as yet, no unstable modes and, 
as far as defects are concerned the field becomes ordered, as before, because 
of the exponential growth of long-wavelength modes, which stop growing 
once the field has sampled the groundstates. What matters is the relative 
weight of these modes in P(k, t) (the 'Bragg' peak) to the fluctuating short 
wavelength modes, since the contribution of these latter is very sensitive to 
the cutoff I at which we look for defects. Only if their contribution to Eq|8| 
is small when field growth stops can a network of vortices be well-defined at 
early times. At t = there is no Bragg peak and the density of line zeroes 
depends entirely on the scale at which we are looking. 



5.3. The Spinodal Time and £ z 



We begin by extending the analysis from t = to later times, still in 
the approximation of a free roll. This needs care for slow quenches since 
the back-reaction serves to hold the field in the vicinity of the intermediate 
ground-states \c/) 2 \ 



6§ (i) where, now 



4>l{t) 



-e(i) 
A 



1 t 



(78) 



for t < tq. 

Prior to the completion of the quench at t 
(]75|), now of the form 



tq, the mode equation 



d 2 , 2 t 

+ k — 

at 1 tq 



xt(t) = 0, 



(79) 



is exactly solvable, as we saw earlier. As in the case of condensed matter 
previously, we have coarse-grained the field by introducing a simple cut-off at 
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k = A = 0(M) in dimensional units, or I = A -1 . In natural units I = O(l). 
The unstable exponentially growing modes appear when 

n 2 k {t) = -k 2 + — > o. (so) 

For fixed k this occurs when t > t k = TQk 2 (or k 2 = k 2 = tjrQ). As before, 
we are not looking for vortices within vortices, but coarse-grain to the cold 
vortex size. 

The WKB solution is adequate for our purposes. G\(r\ t) can be written 
as Gi(r; t) = G exp (r; t) + Gf sc (r; t) where 

G^(r;t)^T fk S *® I 1/3 (S k (t)) + a fe ~ I- 1/3 (S k (t))\ 2 

. | k , , , (81) 

has exponentially growing long wavelength modes and 

Gr(r;t)^T f fk Sk f e* k ->+ J 1/3 (S k (t))-a^ J- 1/3 (S k (t))\ 2 

(82) 

has short wavelength oscillatory modes. In both cases 

S k (t) = f dt' \n k (t')\ = ^_L|i _ t k \ 3 / 2 . (83) 
Jt~ 3^/tq 

We stress that, at the level of a free field (or any self-consistent Gaussian) 
inserting a cut-off is identical to integrating out the short wavelength modes, 
for an unbiased initial state. Provided we are far from the transition we have 
incorporated the initial data into the in (|8l]) and (|82|), which have no A 
or tq dependence. 

For large t the integrand in (^l|) will have its 'Bragg' peak at some 
ko(t) — > as t — > oo, once the angular integrals have been performed. As- 
suming that ko(t) <C h the upper bound in the integral can be dropped 
and \k 2 — k 2 \ approximated by \k 2 \, knowing that there is no singularity at 
k = kf. With nothing to stop \at + a k \ 2 behaving like a nonzero constant 
in the vicinity of k = 0, it can be treated as slowly varying and the integral 
approximated asEl 

G ex P(r;t) oc r^y /2 e ( 4 / 3 )^) 3/2 ^ ^ke^e~ 2 V^ k2 

1 ( _J_\ 3/ \wm/tK) 3 /\-r 2 /2e(t) 



oc ; t——\ e wvwtKr- e -ry^ m (84) 
W)\ 1/2 \^q) 

on performing the k 2 expansion of the exponent, where 

£ 2 W = (85) 
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in natural units tq = 1. Although, like Eq.84, the expression Eq.^ is not 
supposed to be valid for small t, it does embody the Kibble freezeout condi- 
tion in satisfying = 0(1). 

There are two immediate differences between (|84|) and its condensed 



matter counterpart (36). Firstly, in ( p4j) the length £(t) of ( |85|) is both the 
correlation length of the field and exactly the separation of line zeroes ^(t) zero 
as defined through (p6|). Secondly, insofar that the time-dependence of £(t) 
is a guide to early evolution, it grows as i 1 / 4 , rather than the more customary 
i 1 / 2 of condensed matter. 

The calculations above were for a free roll. Let us suppose, provision- 
ally, that the backreaction exerts its influence over such a short time that, 
in effect, it is if it were an instantaneous brake to domain growth. The pro- 
visional freeze-in time t* is then, effectively, the time it takes to reach the 
transient groundstate 4>l(t). That is, G(0;t*) = 0(4>l(t*)), giving 

(v^) 3/2 e~ (4/3)r/fK)3/2 = 0(V 2 (^) 3 / 2 ). (86) 



This gives 

t* _ 1 
Ik 2 



In 



r Q /r 



(1-T G /T C ) 2 , 

and we have used 1—Tg/T c = A, to be specific. For small coupling t*/tjc > 1, 
but will be O(l). A comparison with (|42"|) shows strong similarities with 
condensed matter and that, yet again, t* is insensitive to the parameters of 
the theory. As far as the separation of scales is concerned, we have the same 
effect qualitatively if we had taken t* as the time for the field to reach the 
final ground state as \4> 2 \ = 1/A, rather than the provisional ground states 
(^o(t). Hence we can identify t* with t sp to a good approximation and we 
shall not distinguish between them. 

At this qualitative level the correlation length, and line-zero separation 
length at the spinodal time is 



e(t sP ) = f(t sp ) zer o - e K m 1/2 = o(&)- 



5.4. Fluctuations 

All of this assumes that the oscillatory wavelength terms (|82|) can be 
ignored. Although we have adopted a cutoff at I = O(£o) = 0(M~ 1 ), there 
is still a contribution from modes near M _1 . When we take these modes 
into account the density of line zeroes at t sp can be written as 

Uzero = 72«k(1 + E). (89) 
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The error term E = 0(A 1 / 2 (l/r Q ) 4 / 3 (ln(l/A))- 1 / 3 ) is due to oscillatory 
modes, sensitive to the cut-off A. The condition E 2 1, necessary for the 
line-zero density to be insensitive to scale, is satisfied^ if 

(t q /t ) 2 (1-T g /T c ) <C, (90) 



where C = O(l). This is the QFT counterpart to Eqj6g. If this is the case 



then f 2 = 2ir^Jt sp /tK, that shows yet again that the causal estimate can be 
correct, but for different reasons. 



5.5. Backreaction in QFT 



To improve upon the free-roll result we adopt a mean-field approxima- 
tion along the lines of as we did for the condensed matter systems 
earlier. 

G(r; t) still has the mode decomposition of (76), but the modes \k now 
satisfy the equation 



|_ +k 2 + e(t) + A($2(0))t 



(91) 



where we have taken N = 2 in a large-N (O(N)) theory. 

The end result iscH that, on making a single subtraction at t 



dt 2 



+ k 2 + e(t) + \J d 3 P C(p)lxZ(t) Xp (t) - 1] 



(92) 



which we write as 



5P + k 



e e //(*) 



Xk(t) = 0. 



(93) 



On keeping just the unstable modes in ($ 2 (0))t then, as it grows, its con- 



tribution to (92) weakens the instabilities, so that only longer wavelengths 
become unstable. At t* the instabilities shut off, by definition, and oscilla- 
tory behaviour ensues. 

Explicit calculation shows that the backreaction has little effect for times 
t < t sp . For t > t sp oscillatory modes take over the correlation function and 
we expect oscillations in G(k;t). 

In practice the backreaction rapidly forces e e ff(t) towards zero if the 
coupling is not too small!!. In that case the end result is a new power 
spectrum, obtained by superimposing oscillatory behaviour onto the old 
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spectrum. As a gross oversimplification, the contribution from the earlier 
exponential modes alone can only be to contribute terms something like 



cos k(t — t sp ) + A(k) sin /c(i — t spJ 



2 



(94) 



to G, where the details of A(k) do not concern us. 

The k = mode of Eq.|4| encodes the simple solution Xk=o(t) = a + bt 
when /j, 2 = 0. This has weak causality built into it. Specifically, for r, t — > oo, 

G?(r,i) ~ ~0(2t/r- 1). (95) 
r 

It has to be said that this approximation should not be taken very seriously 
for large t, since we would expect rescattering to take place at times At = 
0(1/ X) in a way that the Gaussian approximation precludes. 

However, it demonstrates that weak causality, implemented by the Gold- 
stone particles of the self-consistent theory, is likely to have little effect on the 
density of line-zeroes that we expect to mature into fully classical vortices, 
depending as it does largely on the behaviour of G at r = 0. 

Finally, all of the above was predicated on the long wavelengths of the 
field having decohered by t sp . It is to this we now turn. 



6. QFT: THE DECOHERENCE TIME 

The empirical TDGL condensed matter Langevin equation ( |30|) that 
we adopted earlier could have been rewritten as a Fokker-Planck equation^ 
for the functional probability density pt [</>]. Our discussion of QFT in the 
previous Section also encoded the corresponding pt [</)]. The difference is 
that, for QFT, we have yet to ascertain whether we can use pt[(j>] to count 
classical configurations, like defects. The reason was given earlier: classical 
defects cannot be identified until adjacent path histories decohere (i.e. until 
quantum interference is negligible) , since only then will pt [(f)] (approximately) 
obey classical laws. 

The mechanism for decoherence is the interaction of the field with its 
environment. In practice, all fields to which the 4>-Reld couple will help it 
decohere. In addition, the short wavelength modes, which play no role in 
field ordering and subsequent defect production, are an important part of 
the environment. For exemplary purposes we consider only the effects of 
short wavelengths of the 0-field itself on decoherence0ii. A more com- 
plete discussion that includes the effects of additional fields has been given 
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elsewhere@. At the level that we have handled QFT so far (a self-consistent 
Gaussian) our cut-off of short wavelength modes is the same as integrat- 
ing them out, and there is no decoherence. We shall need to go beyond a 
free-field approximation to see any effects. 

Decoherence is measured in terms of dissipation. A byproduct of estab- 
lishing decoherence is that we recover a Fokker-Planck-like equation for pt [(/)} 
or, equivalently, a generalised Langevin stochastic equation for the classical 
field. Together with this should be included the explicit dissipation in the 
early universe as a consequence of the expansion of space-time. We shall not 
treat that here, despite its importance (see Refs.39j36,45), 



To show the main ideas we shall consider the simpler system of a single 
real field (whose defects are domain walls) undergoing a very rapid quench. If 
we take the quench time tq ~ tq this is equivalent to taking an instantaneous 
quench, and this we do. There is no difficulty in taking slower quenches in 
principle, it is just that the calculations have yet to be performed. 

That is, yet again we divide the field cp(x) into the long wavelength 
'system'-field </><(x) and the short wavelength 'environment '-field <^>(x), 

73 „, jl/„ -A „«k.x i / \ / j3„, ±!„ -A „ik.x 



</> < (x) = d d X(A(x,i)e lK - x , (j )> (x)= <fa;0(x,i)e lK - x . (96) 

V|k|<l J\k\>l 

= 1 is the dividing line between stable and unstable modes. With tq ~ tq 
) is automatically satisfied and varying the separation wavelength by a 
small amount will have no effect. 

Since it is the system-field 4> < whose behaviour changes dramatically on 
taking T through T c , we adopt an instantaneous quench for T from Tq to 
Tf = at time t = 0, in which e(T) changes sign and magnitude instantly, 
concluding with the value e(t) = —l,t > 0. Meanwhile, for simplicity the 
0> mass is held at the original value eo- 
The classical action separates as 



where 



S[4>] = [</><] + S' cnv [0 > ] + Sint [4><, <f>>], (97) 



4! < l ' 



Smt[<t><, &>] = ~\J d i x(f 1 < {x)(t) 2 > {x). 

We have not included 0<^> and terms in Si nt [4> <1 </>>] since the former 

does not contribute to very long wavelengths and the latter only appears at 
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higher order in A. Our answers will have little to do with the more familiar 
calculations of decoherence, which are motivated by quantum mechanical 
systems with linear coupling to an environment. 

Beginning from an initial thermal distribution peaked around (f) = we 
follow the evolution of the system under the influence of the short wavelength 
environment. The spinodal time is now t sp ~ | ln[l/A]. 

6.1. The Reduced Density Matrix 

Whereas the field probabilities are just the diagonal elements of the 
density matrix p, decoherence is determined from the full density matrix 

0+ 0<, = {^tW <<!>>)■ (98) 

Our Gaussian initial conditions give an uncorrelated thermal initial be- 
haviour p[Tq] = /) < [To]/5>[To] at temperature To. 

The relevant object is the reduced density matrix. It describes the 
evolution of the system under the influence of the environment, and is defined 
by 




The environment will have had the effect of making the system essentially 
classical once p r (t) is, effectively, diagonal. 
Its temporal evolution is given by 




(100) 

where J r is the reduced evolution operator 

ra - ( (101) 

The Feynman-VernoncB influence functional F[(j) + , ] is defined as 

x T >f V<P+ V<p- exp (i{S[<f>$] + S int [cf>+, 0+]) 
J(j)+. J<f> > . 

x exp (-i{S[4>>] +^iat[^<,0>]}), 

where we have used the closed time-path of Fig.l. From the influence func- 
tional we define the influence action 8 A [cf) + ,<fi~] as F ,(/><] = exp id A [0< ,(/><]. 
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We calculate to one loop. After denning 

and using simple and well known identities between propagators, the real 
and imaginary parts of the influence action can be written as 

ReSA = ^- J d A x J d 4 yA(x)K q {x - y)E(y), (102) 

l mS A = ~yJ <^ x J d*yA{x)N q (x,y)A(y), (103) 

where K q (x — y) = lmG^_ + (x,y)6(y° — x°) is the dissipation kernel and 
N q (x — y) = ReG+ + (x,y) is the noise (diffusion) kernel. 



6.2. The Master Equation 

The first step in the evaluation of the evolution equation for p, the 
master equation, is the calculation of the density matrix propagator J r from 



Eq. ( 101 ). We first perform a saddle point approximation 

^r[0< f ,0<f,i|0<i,0<i,io] «expM[c/>+,<y, (104) 

where (f>^ is the solution of the equation of motion |^ + _^_ = with 

boundary conditions <^(io) = <^<i an d cl (t) = <f . Next, we assume that 
the system- field contains only one Fourier mode with k = k^, which we 
identify with the Bragg peak when it arises. 

To estimate tjj it is sufficient to calculate the correction to the usual 
unitary evolution coming from the noise kernel. For clarity we drop the 
suffix / on the final state fields. If A = (</>< 2 — 0< 2 )/2 for the final field 
configurations, then the master equation for p r (t) = (</><|p|<^><) is 

A 2 

i Pr = (4>+\[H, p r ] |</><> - i—VA 2 D{ko, t) Pr + ... (105) 

The volume factor V arises because we are considering a density matrix 
which is a functional of two different field configurations, </><(x) = 0< cos %q.x 
spread over all space. For ko < 1 the diffusion coefficient D(ko,t) shows the 
inevitable exponential growth (To = 0(T C )) 



D(ko,t) ~ T?n(k )exp{2n(k)t}, 



(106) 
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where £l 2 (k) = M 2 — k 2 , (or 1 — k 2 in natural units), associated with the 



instability of the ko mode. Eq. 106 is only valid once 2Q(k)t > 1, which 
will have happened by t = t sp . Once t > t sp the diffusion coefficient stop 
growing, and oscillates around D(ko,t = t sp ). 

On the other hand, for short times i< 1 we find that 



D(k ,t)~T 2 t 



(107) 



up to coefficients 0(1) 



6.3. The Decoherence Time 

We estimate to for the model by considering the approximate solution 



to the master equation (105), 

Pi#<,0<;i] ~ Pr[4><:4><:A ex P 



VT / ds D(ko,s) 
Jo 



(108) 



where p™ is the solution of the unitary part of the master equation (i.e. 
without environment). In terms of the fields 4> = (</>< + </><)/2, and 5 = 
(0+ - r<)l% r = (A 2 /8)^ 2 5 2 . 

The system will decohere when the non-diagonal elements of the reduced 
density matrix are much smaller than the diagonal ones. We therefore look 
at the ratio 



p r [(j) + 5,(p - 5;t] 



p T [<j>,<j)\t\ 



p?$ + 6,$-6;t] 



x exp 



VT I ds D(k ,s) 
Jo 



(109) 



In general, it is not possible to obtain an analytic expression for the ratio 
of density matrices that appears in Eq.(109) since we do not even know the 
diagonal matrix well. However, since diagonalisation will be found to occur 
for t < t sp it is sufficient to neglect the self-coupling of the system field in the 
diagonal matrix elements. In this case the unitary density matrix remains 
Gaussian at all times as 



p^ + 5j-5;t] 



(110) 



where p -1 (£), essentially0@il ((/P)^ 1 , decreases exponentially with time to 
a value 0(A). In the unitary part of the reduced density matrix the non- 
diagonal terms are not suppressed. [This should not be confused with the 
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observation that the unitary Gaussian density matrix does show classical 
correlation, whereby the Wigner functional becomes localised in phase space 
about its classical solutions-^l. However, this has nothing to do with elim- 
inating quantum interference between different field histories.] In order to 
obtain classical behaviour, the relevant part of the reduced density matrix is 
the term proportional to the diffusion coefficient in Eq,( |109| ), since it is this 
that enforces its diagonalisation. 

The decoherence time tr> can be defined as the solution to 

I ps VT f tD ds D(k ,s). (Ill) 
J o 

It corresponds to the time after which we are able to distinguish between 
two different field amplitudes (inside a given volume V). 

To estimate we have to fix the values of V, 5, and <j). V is understood 
as the minimal volume inside which we do not accept coherent superpositions 
of macroscopically distinguishable states for the field. Thus, our choice is 
that this volume factor is 0(M -3 ) = 0(1) in dimensionless units, since the 
Compton wavelength is the smallest scale at which we need to look. Inside 
this volume, we do not discriminate between field amplitudes which differ 
by 0(1), and therefore we take S ~ 0(1). For 4> we set cf> 2 = 0(1/A), its 
post-transition value. 

From the short-times expression for the diffusion coeffient Eq.(107) is 



very easy to show that decoherence does not occur while /jf<l due to the 
small value of the coupling constants. Consequently, in order to evaluate the 
decoherence time in our model, we have to use Eq.( |106| ). We obtain 



exp[2i D ] = 0(1) (112) 

from which it follows that 

1< t D < V (113) 

This result is welcome in that it suggests that our previous analysis can 
probably go ahead to slower quenches. Moreover, if ( |113| ) is valid there we 
anticipate that classical defects can appear as soon as the spinodal time is 
reached. However, some caution is required since perturbation theory is, 
generally, not valid for non-equilibrium processes. Fortunately, within the 
same one-loop framework it is not difficult to find circumstances in which 
flll3l ) holds until t sp . In particular if, in addition to the short wavelength 
modes, we include a large number of weakly coupled external fields to de- 
cohere the long wavelength modes (withoutany linear couplings), the result 
(|113j ) survives. Details are given elsewhereEa. 
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Finally, once the long wavelength modes have decohered they satisy 
stochastic equations of a generalised Langevin 

a 2 0(x) -<j> + ^ 3 (x) + y <j){x) J A K(x, y)4> 2 (y) = \(/>(x) V (x), (114) 

where the Gaussian noise r\ satisfies {r/(x)r)(0)) oc N(x) in this approxima- 
tion. 

Since K is retarded the integral describes Landau damping. Because 
of the instantaneous nature of the quench we shall not pursue this equation 
beyond saying that, even given that, with its weak O(A) multiplicative noise 
and weak quasi-dissipative term, it is very different from the generalised 
relativistic Langevin equation invoked in Refs.|2].||. 



7. EXPERIMENTS: CONDENSED MATTER 

Now that we understand the Kibble-Zurek bounds better, we give a 
brief review of the main experimental evidence, as well as suggesting further 
experiments that could provide further support for them. 

Our first observation is that, in QFT, there is no evidence for or against 
Kibble's bounds. Although we have every confidence in the standard model 
of Electroweak unification (and its breaking) in the early universe, we only 
have very circumstantial evidence for phase transitions at the higher energies 
that are of more interest to us for defect formation. In particular, to date 
there are no reliable sightings of any topological defect that could have been 
produced in the early universe. At best we have possible indirect evidence, 
such as the suggestion that very high energy cosmic rays arise from the in- 
tersection and release of energy from the high energy density cores of cosmic 
strings. [There is a parallel in condensed matter in quasiparticle production 
due to the interaction of vortices in 4 He] . The original hope that large-scale 
structure was determined by cosmic strings (vortices) seems unfounded!!. 

On the other hand, for condensed matter systems there are several ex- 
periments that have been devised to check Zurek's predictions. 



7.1. Superfluid 3 He 

We begin with superfluid 3 He since, to date, it provides the strongest 
support for the validity of Zurek's bounds. 3 He is a fermion but, somewhat 
as in a BCS superconductor, these fermions form the counterpart to Cooper 
pairs. However, whereas the (electron) Cooper pairs in a superconductor 
form a 1 S state, the 3 He pairs form a 3 P state. The order parameter A a i 
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is a complex 3x3 matrix A a i. There are two distinct superfTuid phases, 
depending on how the global SO(3) x £0(3) x U(l) symmetry is broken. 
If the normal fluid is cooled at low pressures, it makes a transition to the 
3 He — B phase, in which A a i takes the form A a i = R a i(uj)e l ® , where R is a 
real rotation matrix, corresponding to a rotation through an arbitrary all. 
Although more complicated than other systems it can be easier to count 
vortices, since one can use NMR to detect them. 

So far, experiments have been of two types. In the Helsinki experiment^ 
superfluid 3 He — B in a rotating cryostat is bombarded by slow neutrons. 
Each neutron entering the chamber releases 760 keV, via the reaction n + 3 
He — > p + 3 He + 760keV. The energy goes into the kinetic energy of the 
proton and triton, and is dissipated by ionisation, heating a region of the 
sample above its transition temperature. The heated region then cools back 
through the transition temperature, creating vortices. Vortices above a crit- 
ical size grow and migrate to the centre of the apparatus, where they are 
counted by an NMR absorption measurement. The quench is very fast, with 
t q/ t o = O(10 3 ). Agreement with Eq|| is good. 

The second type of experiment has been performed at Grenoble and 
LancastercHI. Rather than count individual vortices, the experiment detects 
the total energy going into vortex formation. As before, 3 He is irradiated by 
neutrons. After each absorption the energy released in the form of quasipar- 
ticles is measured, and found to be less than the total 760 keV. This missing 
energy is assumed to have been expended on vortex production. Again, 
agreement with Zurek's prediction Eq.|| is good. 



7.2. 4 He Experiments 

In 4 He the bose superfluid has as its order parameter a complex field 
<f>, whose squared modulus \(j)\ 2 is the superfluid density. The global U(l) 
symmetry breaks to yield simple vortices, whose winding number measures 
the supernow around them. 

The experiments in 4 He, conducted at Lancaster, follow Zurek's origi- 
nal suggestion. A sample of normal fluid helium, in a container with bellows, 
is expanded so that it becomes superfluid at essentially constant tempera- 
ture. That is, we change 1 — T/T c from negative to positive by increasing 
T c . As the system goes into the superfluid phase a network of vortices is 
formed, which are detected by measuring the attenuation of second sound 
within the bellows. A mechanical quench is slow, with tq some tens of mil- 
liseconds, and tq/tq = O(10 10 ). Two experiments have been performed. 
In the firsl£!l fair agreement was found with the prediction Eq.|8[ However, 
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there were potential problems with hydro dynamic effects at the bellows, and 
at the capillary with which the bellows were filled. A second experiment^, 
designed to minimise these and other problems has failed to see any vortices 
whatsoever. 

However, unlike the experiments with s He, the experimental prediction 
for He needs more than causality to go further. The density of vortices n 
is assumed to obey Vinen's equationH 

dn def h 2 

-&r = - X2 m n ^- (115) 

This behaviour, requiring £def{t) °£ t 1 ^ 2 is t ne general behaviour for non- 
conserved order parameters. We stress that causality does not provide a 
value for \2- 

In our simple TDGL model, in which this behaviour follows at early 
times in the absence of thermal fluctuations, X2 = 47r?ir = 4irh/ 'aoTQ. Taking 
To ~ 8.0 x 10~ 12 s and £o ~ the mean-field approximation for 4 He gives 
X2 ~ 5 x Air. [It would have the value 4ir if we motivated To from the Gross- 
Pitaevskii equationlil, in which oqTq = h.] This decay law is assumed in the 
analysis of the Lancaster experiments, in which the empirical value of \2 was 
not taken from quenches, but turbulent flow experiments. It was suggested^ 
that X2 ~ 0.005, orders of magnitude smaller than our suggestion above. 
Although the TDLG theory is not very reliable for 4 He, if our estimate is 
at all sensible it does imply that vortices produced in a temperature quench 
decay much faster than those produced in turbulence. 

As a separate observation, we note that the large value of f 2 in the 
prefactor of 

rizero is, in itself, almost enough to make it impossible to see 
vortices in 4 He experiments, should they be present. 

To compound the problem, for early time at least, thermal fluctua- 
tions are large in the Lancaster experiments. With tq/tq = O(10 10 ) and 
a Ginzburg regime so large that (1 — Tg/T c ) = 0(1) the inequality ( |66| ) is 
hugely violated. In such circumstances the density of zeroes n zero = 0(l~ 2 ) 
after t sp depends exactly on the scale I at which we look and they are not 
candidates for vortices. Since the whole of the quench takes place within the 
Ginzburg regime this is not implausible. However, the Ginzburg-Landau the- 
ory is not reliable for 4 He and, even though the thermal noise never switches 
off, there is probably no more than a postponement of vortex production. 



7.3. Superconductors: Flux Through a Surface 



The situation with local symmetries, in which the 0-field interacts with 
the electromagnetic fields, is different, and only beginning to be understood^. 
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Whereas the 'domains' of the global symmetry breaking are characterised 
by approximately constant field phases, for a local (gauge) symmetry the 
phase can be changed at will by a gauge transformation (except at a zero of 
the field where it is not defined). Phase correlation lengths are not physical 
observables. Nonetheless, vortices are present, identified by their magnetic 
flux. 



For the relativistic case discussed in Ref.53 it is Landau damping that 
controls the relaxation time, commensurate with our earlier comments on 
stochastic equations. Details are given elsewhere in these proceedings!^. In 
this case it is not clear if causality provides a useful constraint, except in 
the case of extremely weak electromagnetic coupling, very fast transitions or 
low temperatures^ when, for condensed matter, we have (0) for the density 
of Abrikosov vorticestJ. 

Even if the situation were clear here, there is a separate issue as to 
how we could confirm any predictions. Since total flux (proportional to the 
number of vortices minus the number of antivortices) is the gauge-invariant 
measurable quantity, rather than the density (the number of vortices plus 
the number of antivortices), we have to infer the separation ^ e f differently. 

We assume a temperature quench of a superconductor with zero aver- 
age flux. Consider a closed path in the superconductor with circumference 
C S> £. The phase difference 9c around the path is gauge invariant and we 
estimate the r.m.s phase difference along the path as 



A6 C = OyC/ttef). (116) 
The variance of magnetic flux within the path is then 




A* c = 0(^/^-1, (117) 



tdef. 

where <E>o = h/2e is the basic unit of magnetic flux. Since topological charge 
is conserved, there is no haste to measure the trapped flux. Such an experi- 
ment has been performed^ with a high-T c superconductor, but there is no 
agreement with (117) from its simple estimates. The reason is unclear. 

In this context we should also mention the production of vortices pro- 
duced in the wakexif Bose-Einstein condensation, for which the Zurek mech- 
anism also applies^, but for which the relevant experiments have yet to be 
performed. 



7.4. Annular experiments 



As an alternative way to check the predictions for ^ef by measuring 
conserved topological charge, Zurek suggested experiments with annuli of 
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condensed matter. If we quenched a superconducting ring of circumference 
C we would expect the variance in flux through it to be given also by Ql 17 ), 



provided its width is is less than £<fe/j and assuming that the fluctuations in 
the electromagnetic field are negligible. 

In fact, it is more convenient to quench annular Josephson Junctions, in 
which two rings of superconductor are held apart by an oxide layer through 
which Cooper pairs can tunnel. Let the complex scalar Higgs fields for the 
superconductors (labeled 1 and 2) be 4>\ and respectively. Despite our 
concerns about the local symmetries within the individual superconductors 
enabling us to gauge transform away the field phases 9\ and 02, the differ- 
ence between these phases determines a physical observable, the Josephson 
current, 

J = J c sm8, (118) 

where 9 = 9\ — 0%. 

The defects of this system are the 'fluxons' of the Josephson Junction 
and are easy to observe experimentally!!. The variance in fluxon number at 
their formation is expected to be 




ANo = — AO » — J^, (119) 



where £z is Zurek's causal length for such a system. There are two competing 
mechanisms. At very early times we can treat the two superconducting am 
nuli as independent. If this decoupling lasts long enough this would suggest^ 
that (A9) 2 = (A#i) 2 + (A#2) 2 , for which we use the Zurek causal constraints 
for free propagators. On the other hand, if the coupling between the JTJs 
is important by the causal time iz, then what matters is that the causal 
horizon can encompass a fluxon. 

Although some caution is required in the interpretation of past experiments 
with JTJs, for which £ ~ C, which were not devised with this prediction in 
mind, we would have expected to see a fluxon a few percent of the time if the 
latter were true. This was the case. If the former mechanism is correct, we 
would have expected to see them more often. New experiments with more 
relevant parameters are being devised and have been reported elsewherelliS. 

In a different context an experiment with an annular superconductor, 
separated into segments by JTJs, shows that the assumed phase separation 
at the onset of a transition (the 'Kibble mechanism' in the context of QFT) 
is correctE3. A further experiment, to quench an annular container of similar 
circumference C of superfluid A He, was also proposed by ZurekS, but has 
not been attempted. Since the phase gradient is directly proportional to the 
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superflow velocity we expect a flow after the quench with r.m.s velocity 

a "=°q/5)- (i2 °» 

Although not large it is measurable, in principle. A similar experiment for 
BEC has been proposed!^, but has also not been performed. 



8. CONCLUSIONS 

The strong causal bounds @ and (0) of Zurek and (|l0|) and (11) of 
Kibble have a common origin that purports to show strong links between 
phase transitions in condensed matter and particle physics. In reality, this 
masks huge differences between the two types of system. Nonetheless, each 
has a certain qualitative validity, but for different reasons. What we have 
found in common is that explicit causality is absent from both. 

In TDGL condensed matter physics we have seen how fluctuations sep- 
arate initially into unstable long wavelength modes which order the field and 
short wavelength thermal fluctuations that make vortices (and other defects) 
fractal. Dimensional analysis sets the Zurek time tz of (||) and length scale 
£,z of (^) as the natural scales when fluctuations can be ignored. However, 
the density of defects is not set by (as the correlation length of the field) 
since it freezes too soon for instabilities to have grown. Instead, the density 
is set by the separation of field zeroes, a very different length, in principle. 

In practice, because instabilities grow exponentially, the time it takes to 
order the field is, qualitatively, 0(tz), with only logarithmic dependence on 
the parameters of the system. Further, since the defect density is, essentially, 
a ratio of moments of the power spectrum of the fluctuations, there is strong 
cancellation of the prefactors that contain the detailed information of the 
system. As a result, if thermal fluctuations are negligible the defect density 
is that of a free field for longer than we might have thought, and Zurek's 
causal predictions are valid, although not for directly causal reasons. If we 
want to see explicit causality in the TDGL formalism we need to go to its 
dual worldline representation^ which we have not discussed here, but which 
we have examined elsewhere^. 

Although thermal fluctuations and long wavelength (unstable) modes 
are not additive in determining defect densities, simple analysis suggests that 
thermal fluctuations are roughly proportional to the current temperature, 
with little memory of the past temperature profile. As a result, the Ginzburg 
regime plays no role in incubating defects. 

Quantum Field Theory has some characteristics in common with con- 
densed matter: exponential (albeit different exponential) behaviour makes 
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t K the natural time-scale, up to logarithmic terms. We recover £x as the 
relevant length-scale after instabilities have formed, rather than before. How- 
ever, there is a potential complication in that we can only adopt a classical 
probabilistic approach to defects once there are no substantial quantum in- 
terference effects. We have shown that decoherence is established before field 
ordering in a simplified model and more work is being done. However, once 
the initial defects are established, qualitatively as suggested by Kibble, their 
long-term behaviour is unclear. Dissipation is determined by Landau damp- 
ing and long wavelength noise is now multiplicative. Unlike the case of the 
TDGL equation, where we have many numerical calculations, we have very 
little for QFT. As before, the Ginzburg regime only looks to make defects 
fuzzy, rather than create them. 

As for prediction and experiment, we have no concrete proposals for 
QFT. For condensed matter, experiments in superfluid 3 -He are in agree- 
ment with Zurek-Kibble, while those for He seem not to be. However, 
we understand this as due, in large part, to an over-optimistic assumption 
about their decay rate. The role of large thermal fluctuations is unclear. Al- 
though the one experiment with High-T c superconductors sees nothing, an 
earlier experiment with Josephson Tunnel Junctions seems to give support 
to Zurek's predictions, and further experiments with JTJs are planned. 

In none of these experiments is there any agreement with the predictions 
of thermal production in the Ginzburg regime. 
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